Enhanced Thermal Stability and Reduced Aggregation in an Antibody Fab Fragment at Elevated Concentrations.

The aggregation of protein therapeutics such as antibodies remains a major challenge in the biopharmaceutical industry. The present study aimed to characterize the impact of the protein concentration on the mechanisms and potential pathways for aggregation, using the antibody Fab fragment A33 as the model protein. Aggregation kinetics were determined for 0.05 to 100 mg/mL Fab A33, at 65 °C. A surprising trend was observed whereby increasing the concentration decreased the relative aggregation rate, ln(v) (% day-1), from 8.5 at 0.05 mg/mL to 4.4 at 100 mg/mL. The absolute aggregation rate (mol L-1 h-1) increased with the concentration following a rate order of approximately 1 up to a concentration of 25 mg/mL. Above this concentration, there was a transition to an apparently negative rate order of -1.1 up to 100 mg/mL. Several potential mechanisms were examined as possible explanations. A greater apparent conformational stability at 100 mg/mL was observed from an increase in the thermal transition midpoint (Tm) by 7-9 °C, relative to those at 1-4 mg/mL. The associated change in unfolding entropy (△Svh) also increased by 14-18% at 25-100 mg/mL, relative to those at 1-4 mg/mL, indicating reduced conformational flexibility in the native ensemble. Addition of Tween or the crowding agents Ficoll and dextran, showed that neither surface adsorption, diffusion limitations nor simple volume crowding affected the aggregation rate. Fitting of kinetic data to a wide range of mechanistic models implied a reversible two-state conformational switch mechanism from aggregation-prone monomers (N*) into non-aggregating native forms (N) at higher concentrations. kD measurements from DLS data also suggested a weak self-attraction while remaining colloidally stable, consistent with macromolecular self-crowding within weakly associated reversible oligomers. Such a model is also consistent with compaction of the native ensemble observed through changes in Tm and △Svh.

Thermodynamic Origin of Differential Excipient-Lysozyme Interactions.

Understanding the intricate interplay of interactions between proteins, excipients, ions and water is important to achieve the effective purification and stable formulation of protein therapeutics. The free energy of lysozyme interacting with two kinds of polyanionic excipients, citrate and tripolyphosphate, together with sodium chloride and TRIS-buffer, are analysed in multiple-walker metadynamics simulations to understand why tripolyphosphate causes lysozyme to precipitate but citrate does not. The resulting multiscale decomposition of energy and entropy components for water, sodium chloride, excipients and lysozyme reveals that lysozyme is more stabilised by the interaction of tripolyphosphate with basic residues. This is accompanied by more sodium ions being released into solution from tripolyphosphate than for citrate, whilst the latter instead has more water molecules released into solution. Even though lysozyme aggregation is not directly probed in this study, these different mechanisms are suspected to drive the cross-linking between lysozyme molecules with vacant basic residues, ultimately leading to precipitation.

Measuring Nonspecific Protein-Protein Interactions by Dynamic Light Scattering.

Dynamic light scattering has become a method of choice for measuring and quantifying weak, nonspecific protein-protein interactions due to its ease of use, minimal sample consumption, and amenability to high-throughput screening via plate readers. A procedure is given on how to prepare protein samples, carry out measurements by commonly used experimental setups including flow through systems, plate readers, and cuvettes, and analyze the correlation functions to obtain diffusion coefficient data. The chapter concludes by a theoretical section that derives and rationalizes the correlation between diffusion coefficient measurements and protein-protein interactions.

Models for Antibody Behavior in Hydrophobic Interaction Chromatography and in Self-Association.

Monoclonal antibodies (mAbs) form an increasingly important sector of the pharmaceutical market, and their behavior in production, processing, and formulation is a key factor in development. With data sets of solution properties for mAbs becoming available, and with amino acid sequences, and structures for many Fabs, it is timely to examine what features correlate with measured data. Here, previously published data for hydrophobic interaction chromatography and the formation of high molecular weight species are studied. Unsurprisingly, aromatic sidechain content of complementarity-determining regions (CDRs), underpins much of the variability in hydrophobic interaction chromatography data. However, this is not reflected in nonpolar solvent accessible surface enrichment at the antigen-combining site, consistent with a view in which hydrophobic interaction strength is dependent on curvature as well as on the extent of an interface. Sequence properties are also superior to surface-based structural properties in correlations with the high molecular weight species data. Combined length of CDRs is the most important factor, which could be an indication of flexibility that facilitates CDR-CDR interactions in mAb self-association. These observations couple to our understanding of protein physicochemical properties, laying the groundwork for improved developability models.

Controlling Phase Separation of Lysozyme with Polyvalent Anions.

The ability of polyvalent anions to influence protein-protein interactions and protein net charge was investigated through solubility and turbidity experiments, determination of osmotic second virial coefficients ( B22), and ζ-potential values for lysozyme solutions.  B22 values showed that all anions reduce protein-protein repulsion between positively charged lysozyme molecules, and those anions with higher net valencies are more effective. The polyvalent anions pyrophosphate and tripolyphosphate were observed to induce protein reentrant condensation, which has been previously observed with negatively charged proteins in the presence of trivalent cations. Reentrant condensation is a phenomenon in which low concentrations of polyvalent ions induce protein precipitation, but further increasing polyvalent ion concentration causes the protein precipitate to resolubilize. Interestingly, citrate does not induce lysozyme reentrant condensation despite having a similar charge, size, and shape to pyrophosphate. We observe qualitative differences in protein behavior when compared against negatively charged proteins in solutions of trivalent cations. The polyphosphate ions induce a much stronger protein-protein attraction, which correlates with the occurrence of a liquid-gel transition that replaces the liquid-liquid transition observed with trivalent cations. The results indicate that solutions of polyphosphate ions provide a model system for exploring the link between the protein-phase diagram and model interaction potentials and also highlight the importance that ion-specific effects can have on protein solubility.

Self-Interaction of Human Serum Albumin: A Formulation Perspective.

In the present study, small-angle X-ray scattering (SAXS) and static light scattering (SLS) have been used to study the solution properties and self-interaction of recombinant human serum albumin (rHSA) molecules in three pharmaceutically relevant buffer systems. Measurements are carried out up to high protein concentrations and as a function of ionic strength by adding sodium chloride to probe the role of electrostatic interactions. The effective structure factors (S eff) as a function of the scattering vector magnitude q have been extracted from the scattering profiles and fit to the solution of the Ornstein-Zernike equation using a screened Yukawa potential to describe the double-layer force. Although only a limited q range is used, accurate fits required including an electrostatic repulsion element in the model at low ionic strength, while only a hard sphere model with a tunable diameter is necessary for fitting to high-ionic-strength data. The fit values of net charge agree with available data from potentiometric titrations. Osmotic compressibility data obtained by extrapolating the SAXS profiles or directly from SLS measurements has been fit to a 10-term virial expansion for hard spheres and an equation of state for hard biaxial ellipsoids. We show that modeling rHSA as an ellipsoid, rather than a sphere, provides a much more accurate fit for the thermodynamic data over the entire concentration range. Osmotic virial coefficient data, derived at low protein concentration, can be used to parameterize the model for predicting the behavior up to concentrations as high as 450 g/L. The findings are especially important for the biopharmaceutical sector, which require approaches for predicting concentrated protein solution behavior using minimal sample consumption.

Coarse-Grained Modeling of Antibodies from Small-Angle Scattering Profiles.

Predicting the concentrated solution behavior for monoclonal antibodies requires developing and using minimal models to describe their shape and interaction potential. Toward this end, the small-angle X-ray scattering (SAXS) profiles for a monoclonal antibody (COE-03) have been measured under solution conditions chosen to produce weak self-association. The experiments are complemented with molecular simulations of a three-bead antibody model with and without interbead attraction. The scattering profile is extracted directly from the molecular simulation to avoid using the decoupling approximation. We examine the ability of the three-bead model to capture features of the scattering profile and the dependence of compressibilty on protein concentration. The three-bead model is able to reproduce generic features of the experimental structure factor as a function of wave vector S(k) including a well-defined shoulder, which is a consequence of the planar structure of the antibody, and a well-defined minimum in S(k) at k ∼ 0.025 Å-1. We also show the decoupling approximation is incapable of accounting for highly anisotropic shapes. The best-fit parameters obtained from matching spherical models to simulated scattering profiles are protein concentration dependent, which limits their applicability for predicting thermodynamic properties. Nevertheless, the experimental compressibility curves can be accurately reproduced by an appropriate parametrization of the Baxter adhesive model, indicating the model provides a semiempirical equation of state for the antibody. The results provide insights into how equations of state can be improved for antibodies by accounting for their anisotropic shapes.

Osmolyte Effects on Monoclonal Antibody Stability and Concentration-Dependent Protein Interactions with Water and Common Osmolytes.

Preferential interactions of proteins with water and osmolytes play a major role in controlling the thermodynamics of protein solutions. While changes in protein stability and shifts in phase behavior are often reported with the addition of osmolytes, the underlying protein interactions with water and/or osmolytes are typically inferred rather than measured directly. In this work, Kirkwood-Buff integrals for protein-water interactions (G12) and protein-osmolyte interactions (G23) were determined as a function of osmolyte concentration from density measurements of antistreptavidin immunoglobulin gamma-1 (AS-IgG1) in ternary aqueous solutions for a set of common neutral osmolytes: sucrose, trehalose, sorbitol, and poly(ethylene glycol) (PEG). For sucrose and PEG solutions, both protein-water and protein-osmolyte interactions depend strongly on osmolyte concentrations (c3). Strikingly, both osmolytes change from being preferentially excluded to preferentially accumulated with increasing c3. In contrast, sorbitol and trehalose solutions do not show large enough preferential interactions to be detected by densimetry. G12 and G23 values are used to estimate the transfer free energy for native AS-IgG1 (Δµ2N) and compared with existing models. AS-IgG1 unfolding via calorimetry shows a linear increase in midpoint temperatures as a function of trehalose, sucrose, and sorbitol concentrations, but the opposite behavior for PEG. Together, the results highlight limitations of existing models and common assumptions regarding the mechanisms of protein stabilization by osmolytes. Finally, PEG preferential interactions destabilize the Fab regions of AS-IgG1 more so than the CH2 or CH3 domains, illustrating preferential interactions can be specific to different protein domains.

Lysine and arginine content of proteins: computational analysis suggests a new tool for solubility design.

Prediction and engineering of protein solubility is an important but imprecise area. While some features are routinely used, such as the avoidance of extensive non-polar surface area, scope remains for benchmarking of sequence and structural features with experimental data. We study properties in the context of experimental solubilities, protein gene expression levels, and families of abundant proteins (serum albumin and myoglobin) and their less abundant paralogues. A common feature that emerges for proteins with elevated solubility and at higher expression and abundance levels is an increased ratio of lysine content to arginine content. We suggest that the same properties of arginine that give rise to its recorded propensity for specific interaction surfaces also lead to favorable interactions at nonspecific contacts, and thus lysine is favored for proteins at relatively high concentration. A survey of protein therapeutics shows that a significant subset possesses a relatively low lysine to arginine ratio, and therefore may not be favored for high protein concentration. We conclude that modulation of lysine and arginine content could prove a useful and relatively simple addition to the toolkit available for engineering protein solubility in biotechnological applications.

Soluble expression of proteins correlates with a lack of positively-charged surface.

Prediction of protein solubility is gaining importance with the growing use of protein molecules as therapeutics, and ongoing requirements for high level expression. We have investigated protein surface features that correlate with insolubility. Non-polar surface patches associate to some degree with insolubility, but this is far exceeded by the association with positively-charged patches. Negatively-charged patches do not separate insoluble/soluble subsets. The separation of soluble and insoluble subsets by positive charge clustering (area under the curve for a ROC plot is 0.85) has a striking parallel with the separation that delineates nucleic acid-binding proteins, although most of the insoluble dataset are not known to bind nucleic acid. Additionally, these basic patches are enriched for arginine, relative to lysine. The results are discussed in the context of expression systems and downstream processing, contributing to a view of protein solubility in which the molecular interactions of charged groups are far from equivalent.

The crystal structure of the FAD/NADPH-binding domain of flavocytochrome P450 BM3.

We report the crystal structure of the FAD/NADPH-binding domain (FAD domain) of the biotechnologically important Bacillus megaterium flavocytochrome P450 BM3, the last domain of the enzyme to be structurally resolved. The structure was solved in both the absence and presence of the ligand NADP(+), identifying important protein interactions with the NADPH 2'-phosphate that helps to dictate specificity for NADPH over NADH, and involving residues Tyr974, Arg966, Lys972 and Ser965. The Trp1046 side chain shields the FAD isoalloxazine ring from NADPH, and motion of this residue is required to enable NADPH-dependent FAD reduction. Multiple binding interactions stabilize the FAD cofactor, including aromatic stacking with the adenine group from the side chains of Tyr860 and Trp854, and several interactions with FAD pyrophosphate oxygens, including bonding to tyrosines 828, 829 and 860. Mutagenesis of C773 and C999 to alanine was required for successful crystallization, with C773A predicted to disfavour intramolecular and intermolecular disulfide bonding. Multiangle laser light scattering analysis showed wild-type FAD domain to be near-exclusively dimeric, with dimer disruption achieved on treatment with the reducing agent dithiothreitol. By contrast, light scattering showed that the C773A/C999A FAD domain was monomeric. The C773A/C999A FAD domain structure confirms that Ala773 is surface exposed and in close proximity to Cys810, with this region of the enzyme's connecting domain (that links the FAD domain to the FMN-binding domain in P450 BM3) located at a crystal contact interface between FAD domains. The FAD domain crystal structure enables molecular modelling of its interactions with its cognate FMN (flavodoxin-like) domain within the BM3 reductase module.

Flavocytochrome P450 BM3 mutant W1046A is a NADH-dependent fatty acid hydroxylase: implications for the mechanism of electron transfer in the P450 BM3 dimer.

Bacillus megaterium P450 BM3 (BM3) is a P450/P450 reductase fusion enzyme, where the dimer is considered the active form in NADPH-dependent fatty acid hydroxylation. The BM3 W1046A mutant was generated, removing an aromatic "shield" from its FAD isoalloxazine ring. W1046A BM3 is a catalytically active NADH-dependent lauric acid hydroxylase, with product formation slightly superior to the NADPH-driven enzyme. The W1046A BM3 K(m) for NADH is 20-fold lower than wild-type BM3, and catalytic efficiency of W1046A BM3 with NADH and NADPH are similar in lauric acid oxidation. Wild-type BM3 also catalyzes NADH-dependent lauric acid hydroxylation, but less efficiently than W1046A BM3. A hypothesis that W1046A BM3 is inactive [15] helped underpin a model of electron transfer from FAD in one BM3 monomer to FMN in the other in order to drive fatty acid hydroxylation in native BM3. Our data showing W1046A BM3 is a functional fatty acid hydroxylase are consistent instead with a BM3 catalytic model involving electron transfer within a reductase monomer, and from FMN of one monomer to heme of the other [12]. W1046A BM3 is an efficient NADH-utilizing fatty acid hydroxylase with potential biotechnological applications.

Acceleration of crystal growth rates: an unexpected effect of tailor-made additives.

The importance of relative growth rates in the preponderance of alpha- over gamma-glycine during solution crystallisation has been confirmed. Most surprisingly tailor-made additives drastically accelerated the growth of gamma-glycine--an unexpected and key factor in the polymorphic outcome of glycine crystallisation.

Density of states for a short overlapping-bead polymer: clues to a mechanism for helix formation?

The densities of states are evaluated for very short chain molecules made up of overlapping monomers, using a model which has previously been shown to produce helical structure. The results of numerical calculations are presented for tetramers and pentamers. We show that these models demonstrate behaviors relevant to the behaviors seen in longer, helix-forming chains, particularly "magic numbers" of the overlap parameter, where the derivatives of the densities of states change discontinuously, and a region of bimodal energy probability distributions, reminiscent of a first-order phase transition in a bulk system.

Structure and aggregation of a helix-forming polymer.

We have studied the competition between helix formation and aggregation for a simple polymer model. We present simulation results for a system of two such polymers, examining the potential of mean force, the balance between intermolecular and intramolecular interactions, and the promotion or disruption of secondary structure brought on by the proximity of the two molecules. In particular, we demonstrate that proximity between two such molecules can stabilize secondary structure. However, for this model, observed secondary structure is not stable enough to prevent collapse of the system into an unstructured globule.